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Abstract. We prove a new lower bound on the indirect Coulomb energy in two dimensional 
quantum mechanics in terms of the single particle density of the system. The new universal 
lower bound is an alternative to the Lieb-Solovej-Yngvason bound with a smaller constant, 
C = (4/3) 3 /V5tt- 1 w 5.90 < C LSY = 192^2^ « 481.27, which also involves an additive 
gradient energy term of the single particle density. 

1. Introduction 

Since the beginning of Quantum Mechanics there has been a wide interest in estimating 
various energy terms of a system of electrons in terms of the single particle density p^{x). 
Given that the expectation value of the Coulomb attraction of the electrons by the nuclei 
can be expressed in closed form in terms of p^,(x), the interest focuses on estimating the 
expectation value of the kinetic energy of the system of electrons and on the expectation 
value of the Coulomb repulsion between the electrons. Here, we will be concerned with the 
latest. The most natural approximation to the expectation value of the Coulomb repulsion 
between the electrons is given by 

1 f 1 

2 J \x-y\ 

which is usually called the direct term. The remainder, i.e., the difference between the ex- 
pectation value of the electronic repulsion and D(p,p), say E, is called the indirect term. In 
1930, Dirac |J gave the first approximation to the indirect Coulomb energy in terms of the 
single particle density. Using an argument with plane waves, he approximated E by 

E « -c D e 2 / 3 / p 4 / 3 dx, (2) 



where cd = (3/4)(3/7r) 1 / 3 rs 0.7386 (see, e.g., [fHJ, p. 299). Here e denotes the absolute value 
of the charge of the electron. The first rigorous lower bound for E was obtained by E.H. 
Lieb in 1979 |]r2fl , using the Hardy-Littlewood Maximal Function ||2"5| . There he found that, 
E > -8.52e 2/3 / p 4/3 dx. The constant 8.52 was substantially improved by E.H. Lieb and S. 
Oxford in 1981 Jl3|, who proved the bound 

E>-Ce 2 / 3 [ p 4 / 3 dx, (3) 



with C = clo — 1-68. The best value for C is unknown, but Lieb and Oxford [13| proved 
that it is larger or equal than 1.234. The Lieb-Oxford value was later improved to 1.636 by 
Chan and Handy, in 1999 0. It is this last constant, as far as we know, that is the smallest 
value for C that has been found to this day During the last thirty years, after the work of 



Lieb and Oxford |13[ , there has been a special interest in quantum chemistry in constructing 

l 
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corrections to the Lieb-Oxford term involving the gradient of the single particle density. This 
interest arises with the expectation that states with a relatively small kinetic energy have a 
smaller indirect part (see, e.g., [1C, and references therein). Recently, Benguria, Bley, 

and Loss obtained an alternative to (^), which has a lower constant (close to 1.45) to the 
expense of adding a gradient term (see Theorem 1.1 in [Q]). 

After the work of Lieb and Oxford ]T3| many people have considered bounds on the indirect 
Coulomb energy in lower dimensions (in particular see, e.g., || for the one dimensional case, 
|16| , [2j| , [^1 and \24\ for the two dimensional case, which is important for the study of 
quantum dots). In this manuscript we give an alternative to the Lieb-Solovej-Yngvason 



bound ||16|| , with a constant much closer to the numerical values proposed in p4]| (see also the 
references therein) to the expense of adding a gradient term. In some sense, the result proven 
here is the analog of the three dimensional result proven in || for two dimensional systems. 

Our main result is the following theorem: 

Theorem 1.1 (Estimate on the indirect Coulomb energy in two dimensions). Let ip G L 2 (R 2N ) 
be normalized to one and symmetric (or antisymmetric) in all its variables. Define 



pip(x) = N I (x, x 2 ,..., x N ) dx 2 ... dx N . 

(JV-l) 



Then, for all e > 0, 

N 



Ety) = ty,£\x i -x j \- 1 il>)-Dfa,Pi,)>-(l + e)P [ pfdx-^l |V P ; /4 | 2 dx 



(4) 



with 



3/2 



y/5n - 1 ~ 5.9045. 



Remarks, i) Our constant /3 ~ 5.9045 is substantially lower than the constant Clsy — 481.27 
found in |16[ (see equation (5.24) °f lemma 5.3 in [[16]^, which is the best bound to date. 



ii) The constant (3 is close to the numerical values (i.e., ~ 1.95 ^ of |23[ (and references 
therein), but is not sharp. 

Our proof relies on a stability result for an auxiliary molecular quantum system in two 
dimensions (which is proven in Section 2) and an observation of Lieb and Thirring |17| . The 
proof of the main theorem is given in Section 3. 

2. A STABILITY RESULT FOR AN AUXILIARY MOLECULAR SYSTEM IN TWO DIMENSIONS 

A key role in our proof of the Lieb-Oxford type bound in two dimensions will be played 
by a stability result on an auxiliary molecular system in two dimensions. This molecular 
system may be viewed as the two dimensional version of the zero mass limit of the relativistic 
Thomas-Fermi-Weizsacker energy functional studied in H (which corresponds to the zero 
mass limit of the model introduced in [|, |7|, §]; the stability properties of the corresponding 
atomic system were studied also in ||). Thus, let us consider the energy functional 



£(p)=a 2 f (Vp 1/A ) 2 dx + b 2 [ p^ 2 dx- [ V(x)p(x)dx + D(p,p) + U, 
Jr 2 Jr 2 ./r 2 



(5) 
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where the potential V is given by 



i=l 1 

which may be viewed as a Coulomb-like potential generated by K point particles (nuclei) of 
(equal) charge z > 0, located at Ri G M 2 (with i = 1,. . .K). Here, the function p(x) > is 
the electronic density of a system of N electrons, and 

If 1 

£) (p ; p) = o/ p( x )] :p(y)dxdy (7) 

1 Jr 2 xr 2 \ x — y\ 

is the electronic repulsion energy. Finally, 

l<i<J<iT 14 J' 

The powers of the first two terms in (|5|), i.e., 

T(p) = a 2 [ (Vp 1/4 ) 2 dx + 6 2 / p 3 / 2 dx, (9) 
Jr 2 Jr 2 

are such that T(p a ) = aT(p), where p a (x) = a 2 p(ax) (with a > 0) is such that L a p a dx = 
J R2 p(x) da;. In other words, the kinetic energy of the electrons scales like one over a length, 
i.e., in the same way as the potential energy. Then, as usual in this situation, the values of 
the coupling constant (i.e., the values of the nuclear charge) will be crucial to ensure stability 
of the system. Our main result in this section is the following stability theorem, which is the 
two dimensional analog of Theorem 1.2 in 

Theorem 2.1. For any a,b > 0, and Ri G M 2 , i = 1,...,K, and for all p > (with 
p G L 3 / 2 (M 2 ) and Vp 1/4 G L 2 (R 2 )), we have that 

i{p) = a 2 [ (V p 1/A ) 2 dx + b 2 [ p 3/2 dx- [ V(x)p(x)dx + D(p,p) + U > 0, (10) 
Jr 2 Jr 2 Jr 2 

where V , D, and U are defined by ffl), (Qj, and respectively, provided, 

0<z<z c (a,b) = ^VT^- (11) 

Here < o < 1 is the only positive root of the quartic equation 

a 2 _ 32(5tt - 1) a 
VT^~ 27 63 M 

on the interval (0, 1). 

In the rest of this section we will give the proof of this theorem, which is similar to the 
proof of Theorem 1.2 in Notice that the upper limit z c (a, b) on z to insure stability is not 
sharp; in other words, there could still be values of z above our z c for which £(p) > 0. We 
start with an appropriate Coulomb uncertainty principle. 

Theorem 2.2. For any smooth function f on the closed disk Dr, of radius R, and for all 

a, b el, we have 

a 2 f \Vf(x)\ 2 dx + b 2 [ f(xfdx>abl ( - - ) f(x) 4 dx. 
Jd r Jd r Jd r \2\x\ RJ 
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To prove the theorem one only imitates the proof of Theorem 2.1 in that deals with 
the three-dimensional case. We start with the following preliminary result which may be of 
independent interest. 



Lemma 2.3. Let u = u(\x\) be a sufficiently smooth real function (so that all the terms in (\Li 
are finite) defined on the interval [0, R], such that u(R) = 0. Then the following uncertainty 
principle holds 



[2u(\x\) + \x\u'(\x\)}f(x) 4 dx 



D, 



<4^jf |V/(z)| 2 dxy /2 ( 



u(\x\) 2 \x\ 2 f(xfdx 



D r 



In ( \13j) there is equality if and only if 



\J^su(s)ds + C' 



1/2 

(13) 
(14) 



for some constants C and A. 

Proof. Set gj(x) = u(\x\)xj. Then we have, 



[2u(\x\) + \x\u'(\x\)]f(x) 4 dx = J2 / ^IM) f{.r)\lr 



Dr 



T I f(x)d J [g j (x)f(x) 3 }dx-3j2 [ f{xfg ] {x)d ] f{x)dx 

i JD R ■ JD R 



= -4 / (Vf(x) } x)u(\x\)f(xYdx. 

JD R 

In the last equality we integrated by parts and made use of the fact that u vanishes on the 
boundary 8Dr. Next, the Schwarz inequality implies 



[2u{\x\) + \x\u'(\x\)}f(x) 4 dx 



Dr 



< 4 



\Vf(x)\ 2 dx 



Dr 



1/2 



u(\x\Y\x\ 2 f(x) b dx 



Dr 



1/2 



In the last expression, equality is obtained if and only if 

djf(x) = ~Xju{\x\)f{xf, 
which after an integration yields the function given by (|i"4|) above. 

Proof of Theorem |^7I . Choosing u(r) = r _1 — R~ l in (p!3|) , we conclude that 



□ 



ab 
~2 



S D (H~£) /(x)4dx - 2ab {L |v/(x)|2dx 



1/2 



1/2 

f{xfdx ) < 



D, 



<a 2 \Vf(x)\ 2 dx + b 2 f(xfdx. 



D f 



D F 



□ 

To prove our main result of this section, i.e., Theorem 2T , we will also need the following 
auxiliary lemma. 
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Lemma 2.4. Let Dl(xo) = {x G M 2 | \x — xo\ < L] and H be a half plane such that 
dist(x , OH) = L and x G H . Then 



\X — X \ 



H\D L (x ) 



Proof. Let us shift the origin of the coordinates to Xq and choose the x Cartesian axes parallel 
to OH. Then in the respective polar coordinates (g, (p), 



— 1 /•7T/2 POO 1 

COS ip \ III 



Rdx = 2/ / — dgdcp + 2 / — dgdip 

IX-Xor J-K/2JL Q Jo Jl Q 

H\D L (x ) 

from which the assertion of the lemma follows by a straightforward integration. □ 

In the sequel we need some notation. We introduce the nearest neighbor, or Voronoi, cells 
|27| (see also the review |T4]|), {Tj}^ =1 , defined by 

Tj = {x \\x - Rj\ <\x - R k \}. (15) 
The boundary of Tj, dTj, consists of a finite number of lines. We also define the distance 

Dj = dist( J R j , dTj) = - min^ - Rj\ \ k ^ j}. (16) 

Finally, we denote by Bj the disk of radius Dj centered at Rj, j — 1, . . . , K. 

One of the key ingredients we need in the sequel is the two dimensional version, (see e.g., 
fT5||), of an electrostatic inequality of Lieb and Yau [18|, [l!|. Define the piecewise function 
$(x) on IR 2 with the aid of the Voronoi cells mentioned above. In the cell Tj, $(x) equals the 
electrostatic potential generated by all the nuclei except for the nucleus situated in Tj itself, 
i.e., for x G Tj, 



i=l 

i+3 



Then, one has (see, e.g., [fT5"|), 

A' 



D(p, p)- f Q(x)p(x) dx + U> Z -J2^T- (18) 
This follows at once from the standard (three dimensional) Lieb- Yau electrostatic inequality 



l8| , |l9f by taking a Borel measure supported on a two dimensional plane, with density p( 



x 



With the help of the Coulomb uncertainty principle and the two dimensional electrostatic 
inequality (|18D, we are ready to prove the following estimate. 



Lemma 2.5. For any p G L 3 / 2 (R 2 ) such that Vp 1/4 G L 2 (R 2 ); for all b x > 0, and b 2 > such 
that bf + b 2 , = b 2 , and z < ab 2 /2, we have 



K -, r 2 



D j 



8 21b\ 



(2z 3 (tt - 1) + vra 3 ^) 



(19) 
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Proof. Setting f{x) A = p(x), splitting K 2 as the disjoint union of the Voronoi cells Tj, using 
Theorem |2.2| in each disk Bj, and discarding the kinetic energy terms (which are positive) in 
the complements Fj \ Bj we conclude that 

t{p)>b{( p 3/2 dx- [ Vpdx + ab 2 J2[ ( 1 -±-)p(x)dx + D(p,p) + U. (20) 

It is convenient to define the piecewise function W(x) as 

+ 1 ~~n~\ = V ( X ) [ixE T 3 \ B J 



x-R 



W ^=^ ( , ab 2 ., (21) 

$(x) + — if x G Bj, 

~? 

Provided z < ab 2 /2 (which we assume from here on), we can estimate from below the sum of 
the second and third integrals in fl2~0f) in terms of W(x) as follows, 

A' 



j=l JB 3 \^ «j| ^7 ~~ 1 Jrj\B j \ x *U\ 



t^jBAx-Ri\ f^J E 



Mi 

A' 



-/ W(x)p(x)dx + J2[ f^- z )P^- dx > 



iy(x)p(x) dx. 



Thus, we can write 

e(p)>ei(p)+6(p), (22) 

with 

£ x (p) = 6 2 f p 3 / 2 dx - / (W- $){x)p(x) dx and, 

f 2 (p) = p) - / $(x)p(x) dx + U. 
Jr 2 

From the definition of £i(p), it is clear that £i(p) > £i(p)> where p(x) = 4{W(x) — ^(x))\/(9bf), 
where as usual u + = max(«, 0). Hence, 



.r- 

where the last equality follows from the definition fl2~T|) of W. As every Tj is contained in a 
half-plane, we may estimate the first integral above with the help of Lemma TA. This way we 

get 

ZM > ~4u - 1) + ™ 3 bl] ^- (23) 
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The lower bound for ^(p) follows at once from (|i~8|), i.e 



N 



3=1 J 



Putting (0), (0), and ( 24]) together the assertion of the lemma immediately follows. 
We end this section with the proof of Theorem I2TT. 



(24) 
□ 



Proof of Theorem ]2j\. Let M{z) stand for the term inside the square brackets on the right 
side of (]l9f) . With z c and a defined by by fllTD and (|l^) respectively, set p = z/z c , and 
62 = pby/l — a. Hence, b\ = b 2 — b\ = b 2 {l — p 2 + p 2 a). Replacing the expressions of b\ and 
62 in the expression for M(z) we get, 

P 2 



where 



M(z) = ^a 2 b 2 (l-a) 



Hp) 



32 a 
P 



h{p)Vl - <t(5tt - 1) 



(25) 



(1 — p 2 + p 2 o) 2 

Here, both p and a belong to the interval [0,1]. It is simple to see that h(p) is strictly 
increasing in the interval [0,1], thus h(p) < h(l) = 1/cr 2 . Using this last inequality in (p5|) 
together with the definition of a, i.e., equation ([12]), we conclude that 

M(z) > 0, 

for all z < z r - □ 



3. Proof of Theorem 1.1 



In this Section we give the proof of the main result of this paper, namely Theorem [L.l| . We 
use an idea introduced by Lieb and Thirring in 1975 in their proof of the stability of matter 
[17| (see also the review article [11] and the recent monograph Hl4|). 

Proof of Theorem [0|. Consider the inequality (|i~0|), with K = N (where N is the number of 
electrons in our original system), z = 1 (i.e., the charge of the electrons), and Ri = Xj (for all 
i = 1, . . ., N). With this choice, according to ([IT]) , the inequality ([LOP is valid as long as a 
and b satisfy the constraint, 

(26) 



2 < aby/1 - a, 



with er G (0, 1) the solution of 



a 



where 



3/2 



a 



(3 2 a 



V57T - 1 ~ 5.9045. 



Then take any normalized wavefunction ip(xi,X2, ■ ■ ■ , x^), and multiply ( |10|) by \i/j(xi, 
and integrate over all the electronic configurations, i.e., on ¥L 2N . Moreover, take p 
We get at once, 

n r r 

E^^^^x.-x^-D^p^y-b 2 pfdx-a 2 |V P ; /4 | 2 da;, 

Jr 2 Jvi 2 



(27) 



(28) 

.,x N )\ 2 
p^(x). 



(29) 
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provided a and b satisfy ( p6| ) and ( p7|) above. Thinking of a G (0, 1) as a free parameter, and 
a, 6 satisfying ( p6| ) and ([27]), and writing £ = (1 — a) /a we get at once from ( |26| ) and fl2~T|) that 

6 2 >(l + e)/3, 

for any e > 0. The theorem then follows by choosing the minimum value of b 2 , i.e., b 2 = (l+e)/3, 
hence a 2 = 4/(/3e). □ 

Remark 3.1. In general the two integral terms in (fjj are not comparable. If one takes a 
very rugged p, normalized to N, the gradient term may be very large while the other term can 
remain small. However, if one takes a smooth p, the gradient term can be very small as we 
illustrate in the example below. Let us denote 

L(p) = [ p(xf /2 dx 
Jr 2 

and 

G(p)= I (Vp{xY^) 2 dx. 
Jr 2 

We will evaluate them for the normal distribution 

p(\ x \) = Ce~ A ^ 2 
where C, A > 0. Some straightforward integration yields 

3 A' 

With C = NA/ir, 

I p(\x\)dx = N, 
Jr 2 

and we have 

G _ 3tt 
~L ~ 2iV' 

i.e., in the "large number of particles" limit, the G term becomes negligible. 
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